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Abstract. Here we investigate the effect lattice geometry has on the lifetime of 
two-dimensional topological quantum memories. Initially, we introduce various lattice 
patterns and show how the error-tolerance against bit-flips and phase-flips depends on 
the structure of the underlying lattice. Subsequently, we investigate the dependence of 
the lifetime of the quantum memory on the structure of the underlying lattice when it 
is subject to a finite temperature. Importantly, we provide a simple effective formula 
for the lifetime of the memory in terms of the average degree of the lattice. Finally, we 
propose optimal geometries for the Josephson junction implementation of topological 
quantum memories. 
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1. Introduction 

Topologically ordered systems are promising candidates for quantum memories. The 
main idea of Kitaev in pQ was to use spin (or qubit) Hamiltonians with topologically 
ordered ground-states to encode quantum information [2j. Logical qubits encoded into 
the ground-state of such systems become virtually resilient against environmentally 
induced local static perturbations to their Hamiltonian. The simplest such memory, 
known as the toric code, is considered as the testbed for studying the properties of 
topological order against various models of errors. For concreteness, we shall restrict 
our study to this model. 

In the presence of probabilistic errors, robustness against perturbations is no longer 
sufficient to ensure self-correction [3j. Assume independent probabilistic bit-flip (Pauli 
X operator) and phase-flip (Pauli Z operator) errors during the preparation of the 
initial state on the physical qubits of the toric code. Retrieval of encoded information is 
possible only if the relative number of X and Z errors is below a certain threshold. This 
error-tolerance threshold has been calculated for the toric code on a square lattices, to 
be approximately 11% [3], as well as for various other regular lattices [5]. 

Consider now a topological memory coupled to a thermal bath. If the system is 
initially in a pure state then the thermal environment will start introducing errors. After 
a certain time one expects that a sufficiently large number of errors will be accumulated 
causing the memory to fail [6]. This time is known as the memory lifetime r jTj. 

Errors in topological systems, such as the toric code, take the form of topological 
defects known as anyons [2j. A finite temperature can cause anyons to be created, 
propagated or annihilated. A fundamental flaw of two-dimensional topological memories 
that rules out self-correction under probabilistic errors is the lack of energy barriers that 
could suppress diffusion of anyons over large distances. For instance, consider a process 
that involves a creation of anyon pair from the ground state, a transport of one anyon 
along a non-contract ible loop on the torus, and a final annihilation of the pair. This 
process enacts a non-trivial transformation on the subspace of ground states. However, 
it can be implemented by a stream of local errors at a constant energy cost [6j |8]. 

In fact, Alicki et al [9] demonstrated that in one and two spatial dimensions, no 
stabiliser Hamiltonian can serve as a self-correcting quantum memory. They showed 
that the relaxation time towards the equilibrium state is a constant independent of the 
lattice size. The physical mechanism behind this instability is the absence of interactions 
between the anyonic thermal excitations. Indeed, [10] demonstrated self-correction 
properties of two-dimensional topological systems by allowing long-range interactions 
between anyons. 

The main goal of this paper is to investigate the possibility of hindering the diffusion 
of anyons by considering the toric code on general, non-regular lattices. By localising 
anyons to certain regions on the lattice, they are less likely to travel far apart, leading 
to better error correction properties. However, deforming the surface code breaks the 
symmetry between the X and Z error correction properties, creating an asymmetry 
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Figure 1: The toric code defined for a general lattice structure. The plaquette, or face operator 
Bf acts non-trivially on the qubits of the face /. The vertex operator A v acts non-trivially on 
the qubits that neighbour vertex v. 

in the error threshold values and their corresponding lifetimes [5]. This is because 
stabiliser operators are no longer symmetric under the lattice duality transformation (i.e. 
exchange of vertices and faces with each other). As an application, we demonstrate how 
such lattice geometries can be employed to enhance topological systems with intrinsically 
biased couplings due to physical implementation such as Josephson-j unction arrays. 

2. The two-dimensional toric code 

The starting point of our study is Kitaev's two-dimensional toric code [lj. Consider a 
general lattice £(E,V, F), as shown in Fig. [TJ with periodic boundary conditions (a 
torus), where E ) V ', and F are sets of edges, vertices, and faces, respectively. A qubit 
is associated with each edge. The stabiliser operators of the code on the lattice C are 
defined for each vertex v G V and face / G F, respectively, as 

A c v = ®Xi,Bf= (g) Z u 

ieE v jeEf 

where Xi and Zi denote Pauli operators on the z-th and j-th qubits respectively, and 
E v and Ef indicate the sets of edges originating from vertex v and surrounding the face 
/ respectively [5]. We define lattice C to be the primal lattice. Its dual lattice, denoted 
by C{E ) V ) F) ) is defined by exchanging the vertices and faces of C with each other. 
The stabiliser space, where information is stored, is composed of states which satisfy 
= W) an d Bf\t/j) = for all v and /. Violations of these stabilizers on vertices 
or faces are associated with different types of anyons. The so called e anyons reside on 
vertices and m anyons on faces. Alternatively, the error syndromes of e and m anyons 
are associated with the vertices on the primal and dual lattices, respectively. Hence, 
the thermal delusion of the e and m anyons depends on the connectivity of the primal 
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(c) Union (d) Dual Union 



Figure 2: The "reduced square" and "union" lattices with their dual lattices. The dual lattice 
is obtained from the primal by exchanging faces and vertices with each other. 

and dual lattices, respectively. Error correction corresponds to a pairwise annihilation 
of anyons which does not form topologically non-contractable loops. 

The toric code has been so far intensively investigated mainly on regular lattices 
[31 |5j. In the case of the square lattice, which is a regular self-dual lattice C = £, 
assuming independent errors, the toric code can be mapped to a random bond Ising 
model, resulting in symmetric error tolerances to e and m anyons and a threshold value 
of 11% [3j. When the toric code is coupled to a thermal bath, one observes a sharp 
transition in time known as the memory lifetime, r. For square lattice, the lifetime of 
the primal and dual lattices is same and for temperature T/J = 0.3 it is approximately 
equals r « 5.8 [7j. 

Our goal is to study the toric code beyond regular self-dual lattices. By considering 
such lattices, one can break the symmetry between properties of X and Z error- 
correction. For example, consider the toric code on the reduced square lattice, as shown 
in Fig. [2} In this case C 7^ L and hence error-tolerances and lifetimes are no longer 
symmetric. Intuitively, e anyons are more localised because the primal lattice is less 
connected leading to a better anyone pairing. On the other hand, m anyons on the dual 
lattice are more sparse because the dual lattice is more connected allowing anyon pairs 
to travel far apart. Indeed, the error-tolerance thresholds of the primal and dual lattices 
are 13% and 8%, respectively. Furthermore, the lifetimes of the lattices are 6.8 and 4.7, 
respectively. 
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Figure 3: (Colour online.) Decay of the logical Z operator in the non-interacting toric code for 
the square lattice (left) and reduced square lattice (right) for grid sizes L = 16, 32, 64, 128 and 
L = 15, 30, 60, 120, respectively. The transition becomes sharper as the system size increases. 
All curves are ensemble averages over 10 4 runs. We have used T/J = 0.3, 7(0) = 7(2 J) and 
applied error correction at the readout stage. 



3. Error model and simulations 

In order to determine the lifetime of the memory against thermal noise, it is necessary 
to simulate the resulting dynamics. In this section we show how to model the classical 
dynamics of the system. We model the interaction of the system with a thermal 
environment by coupling each spin to a bath which can introduce spin errors in the 
initial state |^ ), the ground state of the system [9j [TT]. In the limit of weak coupling, 
the rate for the probabilities to be in state is described by a coupled rate equations 

i 

Here, p^(t) is the time-dependent probability to find the system in the state = 
Tlke^ Zk\ipo) where is the set of affected spin indices in Similarly, p Zi (ip)(t) is the 
probability to be in the state Z^ip). Finally, 7^ = 7(— 0Ji,ip) and 7^ = 7^,^) are 
the transition rates to arrive or leave the state |^), where uj = — e Zi ^ is the energy 
difference between the states and Zi\ip). 

The rates 7(0;) describe transition probabilities between states with energy 
difference w. A standard expression for 7(0;) can be obtained from a spin-boson model 
and reads 



7(0;) = 2k Ti 



e \w\/w c 



1 - e-P" 

where n n is a constant with units l/energy n setting the time scale, ft = with 
T being the temperature of the bath and ks denoting Boltzmann's constant which we 
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will set to one, and uu c is the cutoff frequency of the bath p~21 [13]. For simplicity, 
we assume a large cut off energy uj c — >> 0. For n — 1, the bath is called "Ohmic", 
whereas for n > 2 it is called "super-Ohmic" . In this work we will only consider the 
former case. The relevant rates for our simulation of the toric code are 7(0) (rate for an 
anyon to hop to a free neighbouring site), 7(— 2 J) (rate to create an anyon pair), and 
7 (2 J) = j(—2J)e 2f3J (rate to annihilate a pair of adjacent anyons, obtained from the 
detailed balance property). An important feature of the Ohmic case is that 7(0) 7^ 0, 
thus direct hopping to neighbouring sites is always allowed. 

It is impossible to solve Eq. ([!]) analytically for meaningful system sizes, because 
the number of states p^ grows exponentially with lattice size. We thus stochastically 
simulate the system and obtain the quantities of interest. A single iteration of the 
simulation at time t consists of (i) record the relevant parameters of the system; (ii) 
calculate the total spin flip rate R = £^7(0;^), where ip is the current state of the 
system; (iii) draw the time it takes for the next spin-flip to occur from an exponential 
distribution with rate R; (iv) calculate individual spin-flip probabilities pi = j(uJi^)/R 
and flip a spin at random accordingly. After some initially specified time has been 
reached, we stop and have obtained a single iteration. 

We also apply an error correction scheme at the read out stage. The error 
correction scheme consists of pairing up all detected anyons and then annihilating them 
by connecting each pair with a string of errors from one partner to another. The pairing 
is usually chosen such that all anyons are annihilated with the smallest total number 
of single-spin operators. This is known as the minimal-weight perfect matching and 
can be found in polynomial time with the help of Blossom algorithm [H]. In order to 
find the true matching with minimal weight, one in principle would need to choose the 
set of edges that include all connections from every anyon to every other. This forms 
a complete graph with number of edges that grows quadratically with n. Hence, the 
overall scaling of the matching algorithm becomes formidable for large n. Instead, as a 
good approximation we perform a Delaunay triangulation using Triangle [15], so that 
searching is performed only on anyons that are close to each other resulting in a number 
of edges that is linear in the number of anyons. 

4. Degree as a measure of connectivity 

There are a great deal of different lattices on which the toric code may be defined with 
many distinct characteristics. However, much of their behaviour with respect to the 
propagation of errors depends on a single property, the connectivity of the lattice [5]. 
The average degree of a lattice is a simple and effective measure for the connectivity. 
Lattices with high connectivity will have a large average degree and vice versa. Consider 
a periodic connected planar lattices with a finite number K such that a unit cell of the 
lattice can have vertices of various degrees d 1? rf 2 , ■ ■ • , dx- The degree of a vertex is the 
number of edges originating from the vertex. The average degree of the lattice is defined 
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n=1 n=2 n=3 

Figure 4: Examples of a family of lattices obtained by replacing each edge of the square lattice 
with two and three sub-edges respectively. In general one replaces each edge of the square 
lattice with n-sub-edges. 



as a weighted average 

K 

q = J2 w i d ^ ( 3 ) 

i=i 

where the weights w = {w\, . . . ,Wk) are the proportions in which each node appears 
within a unit cell. Furthermore, by using Euler's formula for planar graphs, one can 
show that for any periodic lattice 

2 + ? = l (4) 
q q 

where q is the average degree of the dual lattice C. For example, the triangular lattice 
has an average degree q — 6, its dual, the hexagonal lattice has an average degree 
q — 3. Intuitively, one expects a highly connected lattice (large q) to be fragile against 
Z errors, resulting in a relatively small error-tolerance threshold. On the other hand, 
by Eq. Q, its dual lattice will be less connected (small q) and more robust against X 
errors, resulting in a relatively bigger error-tolerance threshold. 

To gain an insight into how the threshold probability varies with average degree, 
especially in the q —> oc limit, we will consider the class of lattices given in Fig. [4j 
These are based on the square lattice, but replace each edge with n so-called sub-edges. 
Each edge therefore contains n spins, rather than just one as before, and each vertex 
has a degree of 4n. The effect of a Z error on any of the n spins per edge is identical and 
indistinguishable to that on any other. Each changes the anyonic occupancy of the two 
vertices connected by the edge. A net error on the edge therefore only occurs when an 
odd number of single spin errors are applied, since even numbers of errors cancel each 
other out. If the probability of a single spin Z error is the probability of a net error 
on the edge is therefore 

P(p) = i(l-(l-2p)"). (5) 

Error correction on the n-sub-edge lattice is exactly the same as that on the square 
lattice, except that P is used as the error probability per edge. The threshold probability 
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Lattice 


Average Degree 


Pc T/J= 


0.28 


0.29 


0.3 


0.31 


0.32 


Square 


Primal (4.0) 
Dual (A 0) 


0.1 
1 

U. J. 


7.3 
7 3 


6.4 
fi 4 


5.7 
5 7 


5 

o 


4.6 
4 fi 


Reduced Square 


Primal (3.56) 
Dual (4.57) 


0.13 
0.08 


9.0 
6.1 


7.7 
5.3 


6.8 
4.7 


6.0 
4.1 


5.4 
3.7 


Union 


Primal (4.44) 
Dual (3.63) 


0.09 
0.11 


6.8 

8 


5.8 
7.2 


5.0 
6.5 


4.55 
5.65 


4.15 
5.0 



Table 1: Data for primal and dual, square, reduced square and union lattices. The second 
column lists the average degree of each lattice. The third column lists error-tolerance 
thresholds obtained from simulations. The table also compares the lifetimes of the quantum 
memories as a function of T/J. 



for the n-sub-edge lattice, p™, is therefore that for which P(p™) = p c -> where p c « 0.11 is 
the normal threshold probability for the square lattice. Rearranging Eq. ([5]), we find, 

P n c = \{l - (1 - 2p c ) 1/n ) (6) 
= \{l-{l-2p c )^) (7) 

= I(l_ e |ln(l-2p c)) (g) 

= (-2/g)ln(l-2p c )»0(l/g) (9) 

Here we have used the fact that n — q/A. This shows that the threshold probability will 
decay as 0(1 /q) as q — >> oc. This decay is consistent with that seen in Fig. [5| For the 
dual lattices with q n = 2 + (q — 1/n), a similar argument based on the Chernoff bound 
Po = e - 2 ™(pn-i/2) 2 can ^ e use j to s j low as l/(q — 2) — >► oc, p c « 1/2 — 0(v^ 2). 



5. Numerical results on the relationship between rf, p c and r 

Now that all our tools have been introduced, we may now begin to a study of how the 
error tolerance thresholds and thermal lifetimes of lattices depend on their connectivity. 
We focus on concrete examples of reasonable lattices, i.e. those for which no vertex or 
plaquette operator need act on an nonphysically high number of spins. 

Fig. [5] shows the error-tolerances for various lattices as a function of the degree of 
the lattice. For a square lattice where C = £, the error-tolerances of the primal and 
dual lattices is symmetric. One can break this by considering non-self-dual lattices. 
From Eq. Q, we know that for such lattices, there is an asymmetry in the average 
degrees of the primal and dual lattices. This in turn results in an asymmetry in the 
error-tolerance thresholds of the primal and dual lattices. The data in Fig. [5] confirm 
our earlier intuition that lattices with small average degrees (i.e. less connectivity) will 
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Figure 5: (Colour online.) Left: Simulation data of error-tolerance probability thresholds 
for lattices with various average degrees. We have focused on realistic lattice geometries with 
average degree 3 < q < 6. Right: Analytical data obtained from Eq. ( fl3] ) versus simulation 
data for the lifetime of the quantum memory for various values of T/J and for various lattice 



patterns. This comparison reveals that the Eq. ( 13) provides a good description of the system's 
lifetime. 



have better localisation properties and consequently high error-tolerance thresholds. On 
the other hand, lattices with large average degrees (i.e. low connectivity) will have less 
localisation properties and consequently small error-tolerance thresholds. 

The data in Table (1) show the lifetime of the memory for several lattices at 
different values of T/J. Clearly, the higher the temperature T, the smaller the lifetime. 
Furthermore, for any lattice with C ^ £, the lifetimes of the C and C are not the same. 
Consequently any formula for the lifetime must take into account the structure of the 
underlying lattice, as well as the coupling constant J and the temperature of the system 
T. This because for any non-self-dual lattice the lifetime must exhibits an asymmetry 
and varies with J and T. Let TV be the number of diffusing anyons present in the system, 
the fraction p of spins affected by Z error after a time t is estimated as, 

etN , N 

P=^T> ( 10 ) 

where e is some constant that may or may not depend on the lattice structure. The 
memory fails when p is larger than some critical value p c . This gives a lifetime r for the 
memory 

r = ^(e J / T + l), (11) 

£ 

where we have replaced the density of anyons by the Fermi-Dirac equilibrium distribution 
N/L 2 = (e J / T + l) -1 . It turns out that the constant e is actually independent of the 
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lattice size and structure and that e ~ 1/2. Therefore 

r ^2p c (e J/T + 1). (12) 

This formula gives an excellent approximation for the memory lifetime when the error- 
tolerance threshold is known. Furthermore, one can rewrite this formula by expressing 
p c it in terms of the average degree q 

^2(^) (eJ/r + i) _ (13) 
q 2 

where \i = 0.0231 and v = 0.111 are the fitting parameters obtained from Fig. [5| Thus, 
by simply knowing the average degree of a lattice at a given T / J, one can estimate the 
lifetime of the quantum memory. Fig. [5] compares the lifetimes of the memory obtained 



using Eq. (12) with those obtained by direct simulation. Despite the approximation 



used in deriving Eq. (12), the lifetimes obtained from the two approaches are very close. 

An important conclusion we draw from the above analysis is that for a given value 
of T/J, one can not hope to improve the lifetimes of both the primal and dual lattice, 
relative to the square lattices, solely by changing the lattice geometry. However, under 
a given temperature T, if one consider physical systems with biased couplings J for the 
vertex and the plaquette interactions, then one can find an optimal lattice structure 
that improves the lifetimes of both the primal and dual lattices relative to the square 
lattice under the same temperature T. 

6. Josephson-junction arrays 

In this section we demonstrate that for systems with biased couplings, it is possible to 
improve the coherence time (the minimum of r p and r^) of the quantum memory relative 
to the square lattice by optimising lattice geometry. We demonstrate this by considering 
a Josephson-junction implementation of the quantum memory [16]. The relations that 
relate the charging and the Josephson energies, Eq and Ej ) to the couplings are given 
by 

J v = E3 /4 4 /4 exp(- (? .a yjEj/E c ), (14) 

Jf = y . (15) 

where the subscripts v and / refer to the vertex and face operators respectively, a = 1.61 
is a constant and q and q are related by Eq. Q [16]. Let x be the ratio x = Ej/ E C) for 
x > 1 the system is said to be in the toric code regime. Using this we can simplify the 
above relations 

J v /E c = x 3/4 exp(-q.a y/x), (16) 

J f /E c = -. (17) 
q 

These formulas relate the couplings to the structure of the lattice. We can substitute 



Eqs. (16), (17) in the lifetime formula to obtain formulas for the lifetimes of the primal 
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Figure 6: The coherence time (minimum of r p and r^) of the quantum memory obtained 
from Eqs. (18) and (19) as a function of the average degree for x = 1.01. From the top, 
each curve represent the coherence time at temperature T = 0.002, 0.003, 0.004. The 
circles are the values of the coherence time obtained from direct simulations. 



and dual lattices in terms of q and x 

(18) 

"•TT^ (19) 

Thus, for a fixed x ) by optimising the degree of the lattice q one hopes to be able to 
improve the coherence time of the memory relative to the square lattice. For example, 
using these formulas, for x = 1.01, we calculated the lifetimes as a function of the degree 
for different temperatures, Fig. [6j It is evident that along a given curve, say T = 0.003, 
one can achieve a better coherence time for a Josephson-j unction quantum memory by 
considering, say, a hexagonal lattice (q = 3) rather than a square lattice (q = 4). 

7. Conclusion 

We have studied the effect of deforming lattice geometry on the error-tolerance and 
lifetime of the quantum memory. We have considered the toric code on non-regular 
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lattices obtained by removing or adding edges to the square lattice. Unlike square 
lattice, these lattices are non-self-dual, hence, the connectivity and the error-tolerance 
properties of the primal and the dual lattices are no longer the same. We have shown that 
the average degree of a lattice is a good measure for its connectivity and by considering 
families of lattices we have shown that the two quant it ates are closely related. Lattices 
with small average degree will be more robust against errors and vice versa. We have 
also studied the lifetime of these lattices and provided a formula for calculating the 
lifetime of the memory at a given temperature by using the average degree. Finally, we 
have shown that for systems with biased couplings such as Josephson -junction arrays, 
it is possible to improve the coherence time of the memory relative to the square lattice 
by deforming the lattice geometry. 
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